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Abstract 

We explore scenarios where the highest energy cosmic rays are produced by 
new particle physics near the grand unification scale. Using detailed numeri- 
cal simulations of extragalactic nucleon, 7-ray, and neutrino propagation, we 
show the existence of an interesting parameter range for which such scenar- 
ios may explain part of the data and are consistent with all observational 
constraints. A combination of proposed observatories for ultra- high energy 
cosmic rays, neutrino telescopes of ^ few kilometer scale, and 7-ray astro- 
physics instruments should be able to test these scenarios. In particular, for 
neutrino masses in the eV range, exclusive neutrino decay modes of super- 
heavy particles can give rise to neutrino fluxes comparable to those predicted 
in models of active galactic nuclei. 

PACS numbers: 98.80.Cq, 98.70.Sa, 98.70.Vc, 95.30.Cq 
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I. INTRODUCTION 



The highest energy cosmic ray (HECR) events observed above 100 EeV (1 EeV= 10^^ 
eV) [0,0 are difficult to explain within conventional models involving first order Fermi 
acceleration of charged particles at astrophysical shocks 0. It is hard to accelerate protons 
and heavy nuclei up to such energies even in the most powerful astrophysical objects 0] 
such as radio galaxies and active galactic nuclei. Also, nucleons above ~ 70 EeV lose energy 
drastically due to photo-pion production on the cosmic microwave background (CMB) — the 
Greisen-Zatsepin-Kuzmin (GZK) effect — which limits the distance to possible sources 
to less than ~ 100 Mpc 0. Heavy nuclei are photodisintegrated in the CMB within a few 
Mpc [|^. There are no obvious astronomical sources within ~ 100 Mpc of the Earth. 

A way around these difficulties is to suppose the HECR are created directly at energies 
comparable to or exceeding the observed ones rather than being accelerated from lower 
energies. In the current versions of such "top-down" (TD) scenarios, predominantly 7- 
rays and neutrinos are initially produced at ultra-high energies (UHEs) by the decay of 
supermassive elementary "X" particles related to some grand unified theory (GUT). Such 
X particles could be released from topological defect relics of phase transitions which might 
have been caused by spontaneous breaking of GUT symmetries in the early Universe [Q. 
TD models of this type are attractive because they predict injection spectra which are 
considerably harder than shock acceleration spectra and, unlike the GZK effect for nucleons, 
there is no threshold effect in the attenuation of UHE 7-rays. 

There has been considerable discussion in the literature whether the 7-ray, nucleon, and 
neutrino fluxes predicted by TD scenarios are consistent with all the relevant observational 
data and constraints at various energies |]9HT3|. The absolute flux levels predicted by TD 
models are in general uncertain. While some (though perhaps not all) processes involving 
cosmic strings seem to yield negligibly low fluxes |jT^, other processes such as those involving 
annihilation of magnetic monopole-antimonopole pairs [^,0, cosmic necklaces |]T7|, and 
possible |[l^ (but currently controversial [|19|) direct emission of X particles from cosmic 
strings |2|j2l| can, for reasonable values of parameters, yield X particles at rates sufficient 
to explain the observed HECR flux. 

In this work, instead of trying to calculate the absolute fluxes in specific TD models, 
we use the strategy to numerically calculate the fluxes of nucleons, 7-rays, and neutrinos, 
"optimally" normalize them to match data and constraints, and discuss the feasibility and 
consequences of a set of most "favorable" ranges of the relevant parameters implied by our 
calculations. 

A major new feature of our calculations is that our "all particle" propagation code in- 
cludes the feed-back effect of neutrino cascading on the electromagnetic and hadronic chan- 
nels in a fully self-consistent manner (see below). In addition, spurred by recent experimental 
indications of a possible small neutrino mass, we include in our calculations the effects of 
a small neutrino mass (~ eV) and the consequent Z-boson resonance in the interaction of 
UHE neutrinos with the thermal neutrinos. 
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II. TOP-DOWN MODELS 



The X particles released from topological defects could be gauge bosons, Higgs bosons, 
superheavy fermions, etc. depending on the specific GUT. These X particles would have 
a mass mx comparable to the symmetry breaking scale and would rapidly decay into lep- 
tons and/or quarks of roughly comparable energy. We will accordingly consider several 
possibilities for the decay products. Prior calculations were restricted to decay into only 
quarks. The quarks interact strongly and hadronize into nucleons (iVs) and pions, the latter 
decaying in turn into 7-rays, electrons, and neutrinos. Given the X particle production 
rate, dnx/dt, the effective injection spectrum of particle species a [a = 7,A^, e^,z/) via the 
hadronic channel can be written as {dnx/dt){2/mx){dNa/dx), where x = 2E/mxi and 
dNa/dx is the relevant fragmentation function (FF). For the total hadronic FF, dN^/dx, we 
use solutions of the QCD evolution equations in modified leading logarithmic approximation 
which provide good fits to accelerator data at LEP energies |^ , as well as recently suggested 
extensions for supersymmetry ||2^ (we abbreviate these cases by "no-SUSY" and "SUSY", 
respectively). The difference in the results for these two choices will be a measure of the 
uncertainty associated with the FF. Furthermore, the nucleon content of the hadrons 
is assumed to be in the range 3 to 10%, and the rest pions distributed equally among the 
three charge states (see, however, Ref. [0]). The standard pion decay spectra then give 
the injection spectra of 7-rays, electrons, and neutrinos. The X particle injection rate is 
assumed to be spatially uniform and in the matter-dominated era can be parametrized as 
dnx/dt oc t~'^~^P IQ, where p depends on the specific defect scenario. In this paper we focus 
on the case p = 1 which is representative of a number of specific TD processes involving 
ordinary cosmic strings p5| , p!8| , |20| , pT| , necklaces and magnetic monopoles |T^. Finally, 



we assume that the X particles are nonrelativistic at decay. 



III. NUMERICAL SIMULATIONS 

The 7-rays and electrons produced by X particle decay initiate electromagnetic (EM) 
cascades on low energy radiation fields such as the CMB. The high energy photons undergo 
electron-positron pair production (PP; 77^ e'e"*"), and at energies below ~ 10^"^ eV they 
interact mainly with the universal infrared and optical (IR/0) backgrounds, while above 
~ 100 EeV they interact mainly with the universal radio background (URB). In the Klein- 
Nishina regime, where the center of mass energy is large compared to the electron mass, one 
of the outgoing particles usually carries most of the initial energy. This "leading" electron 
(positron) in turn can transfer almost all of its energy to a background photon via inverse 
Compton scattering (ICS; e% — ^ e'7). EM cascades are driven by this cycle of PP and 
ICS. The energy degradation of the "leading" particle in this cycle is slow, whereas the 
total number of particles grows exponentially with time. This makes a standard Monte 
Carlo treatment difficult. We have therefore used an implicit numerical scheme to solve 
the relevant kinetic equations. A detailed account of our transport equation approach is in 



Ref. ||26[. We include all EM interactions that influence the 7-ray spectrum in the energy 
range 10^ eV < E < 10^^ eV, namely PP, ICS, triplet pair production (TPP; 67;, ee~e~^), 
and double pair production (77;, e^e^e^e^), as well as synchrotron losses of electrons in 
the large scale extragalactic magnetic field (EGMF). 
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Similarly to photons, UHE neutrinos give rise to neutrino cascades in the primordial 
neutrino background via exchange of W and Z bosons p7| , p8| . Besides the secondary neutri- 
nos which drive the neutrino cascade, the W and Z decay products include charged leptons 
and quarks which in turn feed into the EM and hadronic channels. Neutrino interactions 
become especially significant if the relic neutrinos have masses rrii, in the eV range and thus 
constitute hot dark matter, because the Z boson resonance then occurs at an UHE neutrino 
energy i^res = 4 x 10"^^ {eV / m^) eV. In fact, this has been proposed as a significant source of 



HECRs [p9| , |30| . Motivated by recent experimental evidence for neutrino mass we assumed 
a mass of 1 eV for all three neutrino flavors and implemented the relevant W boson interac- 
tions in the t-channel and the Z boson exchange via t- and s-channel. Hot dark matter is also 
expected to cluster, potentially increasing secondary 7-ray and nucleon production pp| , |50 . 
This influences mostly scenarios where X decays into neutrinos only. We parametrize mas- 
sive neutrino clustering by a length scale l^, and an overdensity f^. Values of li, ~ few Mpc 
and fu ~ 20 are conceivable on the local Supercluster scale . 



The relevant nucleon interactions implemented are pair production by protons (p7b 
pe~e^), photoproduction of single or multiple pions (A^7;, N mc, n > 1), and neutron de- 
cay. Production of secondary 7-rays, electrons, and neutrinos by pion decay is also included, 
but is in general negligible in the context of TD scenarios where injection is dominated by 
7-rays and neutrinos over nucleons. We assume a flat Universe with no cosmological con- 
stant, and a Hubble constant oi h = 0.65 in units of 100 km sec~^Mpc~^ throughout. An 
important difference with respect to past work is that we follow all produced particles in the 
EM, hadronic, and neutrino channel, whereas the often-used continuous energy loss (CEL) 
approximation (e.g., [^]) follows only the leading cascade particles. We flnd that the CEL 



approximation can signiflcantly underestimate the cascade flux at lower energies. 

The two major uncertainties in the particle transport are the intensity and spectrum of 
the URB for which there exists only an estimate above a few MHz frequency and the 
average value of the EGMF. To bracket these uncertainties we performed simulations for 
the observational URB estimate from Ref. that has a low-frequency cutoff at 2 MHz 
("minimal"), and the medium and maximal theoretical estimates from Ref. [Q, as well as 



for EGMFs between zero and 10^^ G, the latter motivated by limits from Faraday rotation 
measurements |Q . A strong URB tends to suppress the UHE 7-ray flux by direct absorption 
whereas a strong EGMF blocks EM cascading (which otherwise develops efficiently especially 
in a low URB) by synchrotron cooling of the electrons. 



IV. PARTICLE FLUXES 

We now present results from our flux calculations for a variety of combinations of URBs, 
EGMFs, FFs, fractions of nucleons created in quark fragmentation, and X particle decay 
modes. Tab. | identifles some of the scenarios that were found capable of explaining HECRs 
at least above 100 EeV, without violating any observational constraints, along with the 
predicted composition of the TD component below and above the GZK cutoff. The spectrum 
was normalized in the best possible way to explain observed HECRs as being due either to 
nucleon or 7-ray primaries. The flux below < 20 EeV is presumably due to conventional 
acceleration and was not flt. We remark that above 100 EeV, the best flts for the viable 
scenarios from Tab. |I|have acceptable likelihood signiflcances (see Ref. for details) and are 



4 



TABLE I. Some viable p = 1 TD scenarios explaining HECRs at least above 100 EeV. 
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viable for eV mass neutrinos if their overdensity f,y over a scale obeys specified condi- 
tion for the high URB and vanishing EGMF; for weaker URB/stronger EGMF condition 
relaxes/becomes more stringent, respectively. 

dominant component of "visible" TD flux below and above GZK cutoff at ~ 70 EeV. 

for EGMF > IQ-i" G. 

for EGMF > IQ-^ G. 



consistent with the integral flux above 300 EeV estimated in Refs. , in contrast to direct 
fits to the observed differential flux at 300 EeV |12] which would lead to an overproduction 
of the integral flux at higher energies. 

Figs. show the fluxes of some scenarios indicated in Tab. |, along with current 
observational constraints on and projected sensitivities of some future experiments to 7- 
ray and neutrino fluxes. This demonstrates consistency with present constraints within the 
normalization ambiguity. In particular, EM energy injected at high redshifts is recycled by 
cascading to lower energies, as can be seen in Fig. ^ TD models are therefore significantly 
constrained [p|,p!0| by current limits on the diffuse 7-ray background between 30 MeV and 100 
GeV 1^ which acts as a "calorimeter" and requires Q^u ~ 2 .2 xlO'"^^ h{3p—l) eVcm"^ sec~i 
for the total energy injected into the EM channel. On the other hand, it is not clear whether 
the observed diffuse background above 10 GeV can be fully accounted for by conventional 
sources such as unresolved blazars [O] and it has been suggested that decays of heavy 
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FIG. 1. Energy spectra of nucleons, 7-rays and neutrinos for the TD model with mx = 10^^ 
GeV, p = 1, and the decay mode X ^ q + q, assuming the high URB version and an EGMF of 
10^^*^ G. Thick and thin lines represent the SUSY and no-SUSY FFs, respectively. 1 sigma error 
bars are the combined data from the Haverah Park ||3^, Fly's Eye Q and AGASA |||] experiments 
above 10 EeV. Also shown are piecewise power law fits to the observed charged cosmic ray flux 
below 10 EeV, the EGRET measurement of the diffuse 7-ray flux between 30 MeV and 100 GeV 
p7| , and experimental neutrino flux limits from Frejus [^] and Fly's Eye |3^, as well as projected 
neutrino sensitivities of the future Pierre Auger [40| and NASA's OWL |42| projects. 



particles may provide a significant contribution in this energy range [20]. As can be seen in 
the figures, this is also the case for the TD scenarios studied here. In these scenarios, the 
CMB depletes the photon flux above 100 TeV, and the IR/0 background in the range 100 
GeV-lOOTeV, recycling it to energies below 100 GeV (see Fig. |l|). The resulting background 
is not very sensitive to the specific IR/0 background model, however Constraints from 
limits on CMB distortions and light element abundances from ^He-photo disintegration are 



comparable to the bound from the directly observed 7-rays |T0 



Figs. and ^ compare the UHE fluxes from four TD scenarios indicated in Tab. |I[ Fig. |^ 
compares the SUSY and no-SUSY FF for nix = 10^^ GeV and is just the UHE part of Fig. |l|. 
Fig. I compares the two decay channels X — > g + / and X ^ q + q ior mx = 10^^ GeV, 
assuming the no-SUSY FF. Both figures assume the high URB, an EGMF of 10^^'^ G and 
^ 10^^^ G, respectively, and a fraction fj^ — 10% of nucleons created in quark fragmentation. 
The present energy injection rate Qhecr required to produce the UHE fluxes ja{E) can be 
estimated as 
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FIG. 2. Blow up of Fig. |l| for the fluxes at energies above 1 EeV. The tau neutrino fluxes were 
omitted for clarity. 
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where x = 2E/mx, \a{E) is the effective attenuation length of species a, and the fiducial 
values are for E = 100 EeV, and the SUSY FF for mx = 10^^ GeV. For the SUSY and 
no-SUSY FF, Qhecr turns out to be minimal around mx ~ 10^^ GeV and 10^^ GeV, respec- 
tively, and increases below and above that. This is confirmed by the numerical calculations, 
as can be seen from Figs. ^ and ^ and from Tab. |. We therefore conclude that for most 
combinations of the URB and the EGMF, the most poorly known astrophysical ingredients, 
one can find combinations of possible decay modes and FFs that make p = 1 TD models 
with homogeneous source distribution viable HECR explanations for 10^^ GeV ^ mx ~ 10^^ 
GeV. We note in this context that in some GUT models, certain baryon number violating 
decay modes involving leptons and quarks may violate limits on proton decay if mx is too 
far below 10^^ GeV, and may therefore be disfavored, see, for example Ref. |45[ . 

The energy loss and absorption lengths for UHE nucleons and photons are short (^100 
Mpc). Thus, their predicted UHE fluxes are independent of cosmological evolution. The 
7-ray flux below ~ 10^^ eV, however, scales as the total X particle energy release integrated 
over all redshifts and increases with decreasing p [|l^ roughly as l/(3p — 1). Scenarios with 
p < 1 are therefore in general ruled out (see Figs. |I| and |^), whereas constant comoving 
injection rates {p = 2) are well within the limits. Since the EM flux above ~ 10^^ eV is 
efficiently recycled to lower energies, the constraint on p is in general less sensitive to mx 
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FIG. 3. Same as Fig. ^, but for mx = 10^'^ GeV, and the no-SUSY FF, assuming a vanishing 
EGMF. Here, the thick and thin hnes represent the decay modes X ^ q + q and X ^ I + q, 
respectively. The same normahzation of the GeV 7-ray flux as in Figs. |2| was used. 



then expected from earlier CEL-based analytical estimates 

A specific p = 2 scenario is realized in the case where the supermassive X particles have 
a lifetime longer than the age of the Universe and constitute part of cold dark matter, for 
which non-thermal production in the early Universe has recently been identified as a serious 
possibility [^. In this case, local clustering of the sources in the galactic halo has to be 



taken into account which provides the dominant contribution to observable fluxes [0. As 
a consequence, predicted spectra and composition just reflect the injection spectrum, and 
the diffuse 7-ray background at EGRET energies is not a serious constraint. 

We now turn to signatures of TD models at UHEs. The full cascade calculations predict 
7-ray fluxes below 100 EeV that are a factor ~ 3 and ~ 10 higher than those obtained 
using the GEL or absorption approximation often used in the literature in the case of 
strong and weak URB, respectively. This is also reflected by comparing Eq. (|l]) for the 7-ray 
flux with the energy injection rate Qem allowed by the EGRET observations, which yields 
A-y ~ 100 Mpc. Again, this shows the importance of non-leading particles in the development 
of unsaturated EM cascades at energies below ~ 10^^ eV. As a consequence, in all viable 
HEGR explaining cases with only quarks among the X particle decay products, we obtain 
7/nucleon ratios above 200 EeV that are ^ O.I/Zat for nix ~ 10^^ GeV, and about a factor 
2 smaller for mx ^ 10^^ GeV, even for the maximal URB, if the EGMF is < 10"" G. This 
ratio is about a factor 3 higher for the decay modes containing a charged lepton. Although a 
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7-ray primary for the HECR events is somewhat disfavored currently , the compositional 
issue is not settled yet, but future experiments such as the Pierre Auger project ||5^ should 
be able to distinguish 7-ray and nucleon primaries and test this signature. We stress that 
there are viable scenarios with nucleon fluxes that are comparable with or even higher than 
the 7-ray flux at all energies in case of the high URB and/or for a strong EGMF, and 
/tv — 10%, see Figs. |^ and and Tab. |. The predictions from the SUSY FF in Fig. ^ even 
seems able to explain all cosmic rays above ~ 50 EeV, including the dip around 100 EeV, 
as a cross-over from nucleon domination to an about equal mixture of 7-rays and nucleons. 
The low rrix, pure quark decay modes such as the one shown in Fig. | may be able to 
explain all cosmic rays above 10 EeV by nucleon primaries, but also tend to produce a more 
rapid fall-off of fluxes beyond 100 EeV, which constitutes another testable signature. The 
7/nucleon ratio above 100 EeV is about a factor 5 and 10 higher in the medium and minimal 
URB, respectively, as compared to the strong URB case, and in general decreases strongly 
with increasing EGMF > 10"^^ G. 
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FIG. 4. Same as Fig. |l|, but for the pure neutrino decay mode with no EGMF. Shown are the 
maximal UHE neutrino fluxes allowed by the EGRET limit for mx = 10^^ GeV (thick lines) and 
nix = 10^^ GeV (thin lines). For neutrino clustering the lower limits from Tab. |, required to 
explain HECRs, were assumed. This would correspond to overdensities of ~ 30 and ~ 75 over a 
scale lu — 5 Mpc. 



As indicated in Tab. |, another interesting scenario involves the pure neutrino decay 
modes, also shown in Fig. ^for mx < 10^^ GeV. Here, the 7-rays and nucleons are produced 
as secondaries from the interactions of these UHE neutrinos with the relic neutrinos. Because 
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7-rays and nucleons above 100 EeV must have been produced within a distance Xa — few 
Mpc from the observer, their flux is dominantly produced by interactions with the locally 
clustered neutrinos if li,fu ^ Xa- In this case, the energy fluence in the secondaries is about 
fz{fJu/Xz) times the energy fluence in primary neutrinos around the Z resonance, where 
fz ~ 3% is the fractional width of the Z and Xz — 38 Gpc is the neutrino mean free path at 
the Z resonance at zero redshift. In contrast, at energies where the Universe is transparent 
for particles today, the dominant contribution to their production by UHE neutrinos comes 
from interactions with the unclustered relic neutrino component at high redshift. This is 
because for energies E ^ -Eres, the probability for both resonant and non-resonant interaction 
with the relic background per redshift interval is roughly (1 + zY^'^fzto/Xz in the matter 
dominated regime, where to is the age of the Universe (for E < E^^s this probability decreases 
linearly with E). Because the Universe is opaque for 7-rays above ~ 100 TeV, this implies 
that the diffuse 7-ray background below some energy E is sensitive to the injection history 
at z ^ (100 TeVZ-E)^/^. This explains why the 7-ray background is steeper below 10 GeV 
than in the scenarios where its dominant production is not by neutrino interactions, see 
Fig. ^. As a further consequence of neutrino interactions, the secondary neutrino fluxes 
below an energy E < i^res are sensitive to the injection history at z ^ E^-cs/E. For p = 1 
scenarios, all other fluxes are insensitive to the injection history at z ^ 100. Since we are 
mainly interested in neutrino fluxes above 10 EeV and 7-ray fluxes above 100 MeV, it was 
therefore sufficient to integrate injection up to z = 10^ which also approximately marks the 
transition to radiation domination. In addition, for p = 1, the scaling of neutrino interaction 
rates implies that the energy content in the secondaries, and thus in particular in the low 
energy cascade 7-rays, constitutes a few percent of the energy in UHE neutrinos. This fixes 
the maximally allowed UHE neutrino flux which is shown in Fig. ^ and implies the lower 
limit on l,yfi, given in Tab. | which is required if secondaries of UHE neutrino interactions 
are to explain HECRs. The maximal energy injection rate in neutrinos today allowed by the 
EGRET limit is correspondingly higher than the upper bound on Qem by about a factor 
100. Observational consequences of the UHE neutrino fluxes are discussed in the following 
section. 

The spectra predicted by scenarios where the X particles decay into more than two 
quanta are qualitatively similar to the ones for decay into two particles of the same type. 
The details, however, depend on the energy distributions of the decay products. To avoid 
introducing further model dependent parameters, we do not consider such refinements in the 
present paper as we do not consider scenarios where the X particles themselves are created 
with relativistic energies. 

V. NEUTRINO FLUX DETECTION 

In order to discuss the prospects of detectability of neutrino fluxes in TD scenarios we 
express the (in general energy dependent) experimental sensitivities in terms of the ice or 
water equivalent acceptance A{E) (in units of volume times solid angle). Future neutrino 
telescopes of kilometer scale or larger will utilize the detection of Cherenkov radiation from 
muons and EM showers created in interactions of the neutrinos with nucleons either in 
ice or in the deep sea. Examples for experiments that aim at this effective size are the 
ICECUBE version of the AMANDA experiment at the South Pole, as well as the Radio Ice 
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Cherenkov Experiment (RICE) that aims at measuring the radio pulse from the neutrino 
interaction, the French Astronomy with a Neutrino Telescope and Abyss environmental 
RESearch (ANTARES) proposal, and the NESTOR project in the Mediterranian. 

An alternative method is to search for extensive air showers initiated by electrons pro- 
duced by neutrinos via the charged current process. The interaction length of cosmic ray 
hadrons and 7-rays is ~ 100 g cm~^ above 10 EeV and the probability of these strongly in- 
teracting particles initiating air showers deeper than 1500 g cm~^ is negligibly small. Thus, 
showers starting deep in the atmosphere must be produced by penetrating particles such 
as neutrinos. Large neutrino detectors for measuring HECR air showers using the air fluo- 
rescence technique, such as the High Resolution Fly's Eye now under construction ||5T| or 



the planned Japanese Telescope Array will have the potential to search for deeply pen- 



etrating showers (DPS) initiated by neutrinos p8[. Their resolution of measurement of the 
atmospheric depth at which the shower has its maximum particle density is expected to be 
less than 30 g cm~^ and the discrimination between DPS and the regular air showers would 
be relatively straightforward. A possible contamination by a potential background of DPS, 
secondary showers that result from tau lepton decays deep in the atmosphere or from 7-ray 
bremsstrahlung by muons, has been estimated to be less than 10~^ for 10 years observation 
by a typical fluorescence detector. Hence UHE neutrino astronomy with air fluorescence 
detectors is not background limited p8 |. 



In addition, a giant surface array such as the proposed Pierre Auger project |^ also 
has significant sensitivity for neutrino detection by search for horizontal air showers. The 
recently proposed satellite observatory concept for an Orbiting Wide-angle Light collector 
(OWL) would increase the sensitivity to horizontal air showers by at least another order 
of magnitude. 

Detection rates can be obtained by folding the predicted fluxes with the product of the 
charged current neutrino-nucleon cross section for which we use the recent parametrization 
(T^NiE) ~ 2.82 X 10-32(^/10 EeV)°-^°2^,jn2 ^^^^ ^-^^ acceptance A{E). Since the astro- 
physical "background" from other sources of UHE neutrinos, most notably active galactic 
nuclei and Gamma Ray Bursts , and the secondary neutrinos produced by photopion 

production by HECR [^, is expected to be negligible above 10 EeV, we present integral 



event rates for neutrinos above 10 EeV in Fig |^ for the viable HECR explaining TD models 
from Tab. |. We furthermore assume an acceptance scaling as A{E) oc E^'"^^ which seems to 
be imphed by experimental studies. 

For a given mx, the maximum of the neutrino event rates over all decay modes except the 
ones only involving neutrinos is typically reached for the pure quark decay modes, except for 
mx = 10^3 GeV, where the lu mode produces the highest rates. As can be seen from Fig. |^, 
for all nix this maximum actually saturates the general bound on the integral neutrino 
detection rate R{E) pointed out in Ref. [113], 



R{E) < 0.34; 



2-iT km'^ sr 



77. X -0.6 

yr-S (2) 



10^9 eV 



for E ^ 1 PeV, where r is the ratio of energies injected into the neutrino versus EM channel. 
This is not surprising because for all decay modes except the ones only involving neutrinos, 
r < 0.5. The constraint Eq. (|^) is independent of the FF and arises from comparing the 
energy content in neutrinos and 7-rays, the latter being bounded from above by the EGRET 
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FIG. 5. Maximal event rates for muon neutrinos and anti-neutrinos in a detector of 27rkm'^ sr 
acceptance for the viable scenarios from Tab. ||, ordered by row number for given mx- Electron 
neutrino event rates are about a factor 2 smaller. The rates for tau neutrinos are at least a factor 
100 smaller still, except if produced directly in the decay. The telescope array is roughly sensitive 
to the range above the horizontal line, assuming a duty cycle of 10 % and a lifetime of 10 years. 



measurement. 

The highest possible rates are reached for the exclusive neutrino decay mode at mx = 
10^^ GeV for which the bound Eq. (Q) is not applicable because r = oo, and the relevant 
quantity is the fraction of energy produced as secondary 7-rays instead. As can be seen 
from Fig. ^ the neutrino flux continues down to ~ 10^^ eV in these scenarios and can be 

The maximally 



comparable to fluxes predicted by models of active galactic nuclei ||5^j5^ 
possible event rates from muon neutrinos above 1 PeV per year in a 2ii km^ sr detector are 
~ 5.5 for mx = 10^^ GeV, and ~ 3.5 for mx = 10^^ GeV. 

In general, we conclude that at least the highest rates predicted by TD models should 
be observable by next generation experiments such as the Pierre Auger Observatory and 
especially the OWL project, as can also be seen from the sensitivities shown in the figures. 



VI. SUMMARY 



Apart from the decay spectra and rates, the uncertainty of flux predictions in TD sce- 
narios is governed by astrophysical uncertainties, mainly the universal radio background and 
the large scale extragalactic magnetic field. Our calculations show, however, that for most 
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combinations of likely values for these astrophysical parameters and the energy scale of new 
physics, there are possible decay modes and fragmentation functions that lead to scenarios 
explaining the highest energy cosmic rays above the GZK cutoff, and some of them even 
down to ~ 10 EeV, without violating observational constraints on 7-ray and neutrino fluxes. 
For example, an X particle of mass nix — 10^^ GeV decaying into quarks with a fragmenta- 
tion function motivated by supersymmetry can explain cosmic rays above ~ 50 EeV. This 
scenario predicts a transition from a nucleon dominated component to an about equal mix- 
ture of nucleons and 7-rays above ~ 100 EeV in case of a relatively strong universal radio 
background and a large scale magnetic held ^ 10~^° G, a signature that should be testable 
within the next few years. Other tests involve GeV 7-rays whose flux comes close to the 
EGRET measurement, and ultra high energy neutrino fluxes that should be detectable by 
^ few km scale neutrino observatories which are now in the planning stage. 

Another interesting viable class of scenarios involves pure neutrino decay modes in the 
context of eV neutrino masses which can yield even higher neutrino event rates up to a 
few per year in km scale detectors above ~ 10 EeV for mx ^ 10^^ GeV. The neutrino flux 
extends down to ~ 1 PeV in these models where it can be comparable to predictions from 
models of active galactic nuclei. Furthermore, for a modest amount of clustering of neutrino 
dark matter on the scale of the local Supercluster, secondary 7-ray and nucleon production 
by neutrino interactions with the clustered component can provide a significant fraction of 
the highest energy cosmic ray flux. 

ACKNOWLEDGEMENTS 

Special thanks go to the late David Schramm for his constant encouragement and sup- 
port for interdisciplinary research in particle astrophysics. We also thank Wolfgang Ochs, 
Jim Cronin, Chris Hill, and Felix Aharonian for stimulating discussions, Paolo Coppi for 
collaboration in earlier work, and Haim Goldberg for helpful correspondence. P.B. is sup- 
ported at NASA-GSFC by NAS/NRC and NASA. At the University of Chicago this work 
was supported by DOE, NSF and NASA. 



13 



REFERENCES 



D. J. Bird et al, Phys. Rev. Lett. 71, 3401 (1993); Astrophys. J. 424, 491 (1994); ibid. 
441, 144 (1995). 

N. Hayashida et al, Phys. Rev. Lett. 73, 3491 (1994); S. Yoshida et al, As- 
tropart. Phys. 3, 105 (1995); M. Takeda et al, Phys. Rev. Lett. 81, 1163 (1998). 
for a review see, e.g., R. Blandford and D. Eichler, Phys. Rep. 154, 1 (1987). 
A. M. Hillas, Ann. Rev. Astron. Astrophys. 22, 425 (1984). 

K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. T. Zatsepin and V. A. Kuzmin, Pisma 
Zh. Eksp. Teor. Fiz. 4, 114 (1966) [JETP. Lett. 4, 78 (1966)]. 

G. Sigl, D. N. Schramm, and P. Bhattacharjee, Astropart. Phys. 2, 401 (1994); J. W. El- 
bert and P. Sommers, Astrophys. J. 441, 151 (1995). 

J. L. Puget, F. W. Stecker, and J. H. Bredekamp, Astrophys. J. 205, 638 (1976); 
L. N. Epele and E. Roulet, e-print |astro-ph/9808104| ; F. W. Stecker and M. H. Salamom, 



e-print |astro-ph/9808110 



P. Bhattacharjee, C. T. Hill, and D. N. Schramm, Phys. Rev. Lett. 69, 567 (1992). 
X. Chi et al, Astropart. Phys. 1, 129 (1993); iHd. 1, 239 (1993). 

G. Sigl, K. Jedamzik, D. N. Schramm, and V. Berezinsky, Phys. Rev. D 52, 6682 (1995). 
R. J. Protheroe and P. A. Johnson, Astropart. Phys. 4, 253 (1996). 
R. J. Protheroe and T. Stanev, Phys. Rev. Lett. 77, 3708 (1996); erratum Phys. Rev. 
Lett. 78, 2420 (1997). 

G. Sigl, S. Lee, D. N. Schramm, and P. S. Coppi, Phys. Lett. B 392, 129 (1997). 

A. J. Gill and T. W. B. Kibble, Phys. Rev. D 50, 3660 (1994). 

C. T. Hill, Nucl. Phys. B 224, 469 (1983). 

P. Bhattacharjee and G. Sigl, Phys. Rev. D 51, 4079 (1995). 

V. Berezinsky and A. Vilenkin, Phys. Rev. Lett. 79 (1997) 5202. 

G. Vincent, N.D. Antunes, and M. Hindmarsh, Phys. Rev. Lett. 80, 2277 (1998). 

J. N. Moore and E. P. S. Shellard, e-print |hep-ph/98085^. 



P. Bhattacharjee, Q. Shafi, and F. W. Stecker, Phys. Rev. Lett. 80, 3698 (1998). 
U. F. Wichoski, J. H. MacGibbon, and R. H. Brandenberger, e-print |hep-ph/9805419 



Yu. L. Dokshitzer, V. A. Khoze, A. H. Miiller, and S. L Troyan, Basics of Perturbative 
QCD (Editions Frontieres, Singapore, 1991). 

V. Berezinsky and M. Kachelriess, e-print |hep-ph/9803500| , to appear in Phys. Lett. B. 



M. Birkel and S. Sarkar, e-print |hep-ph/ 9804285 



P. Bhattacharjee and N. C. Rana, Phys. Lett. B 246, 365 (1990). 
S. Lee, Phys. Rev. D 58, 043004 (1998). 

T. J. Weiler, Phys. Rev. Lett. 49, 234 (1982); Astrophys. J. 285, 495 (1984). 

S. Yoshida, Astropart. Phys. 2, 187 (1994); S. Yoshida, H. Dai, C. C. H. Jui, and 

P. Sommers, Astrophys. J. 479, 547 (1997). 

T. J. Weiler, e-print |hep-ph/ 9 710431] , to appear in Astropart. Phys. 

S. Yoshida, G. Sigl, and S. Lee, e-print [hep-ph / 9808324] , submitted to Phys. Rev. Lett. 



F. A. Aharonian, P. Bhattacharjee, and D. N. Schramm, Phys. Rev. D 46, 4188 (1992). 
T. A. Clark, L. W. Brown, and J. K. Alexander, Nature 228, 847 (1970). 

R. J. Protheroe and P. L. Biermann, Astropart. Phys. 6, 45 (1996). 
P. P. Kronberg, Rep. Prog. Phys. 57, 325 (1994). 

G. Sigl, S. Lee, D. N. Schramm, and P. Bhattacharjee, Science 270, 1977 (1995). 



14 



[36] M. A. Lawrence, R. J. O. Reid, and A. A. Watson, J. Phys. G Nucl. Part. Phys 17, 733 
(1991). 

P. Sreekumar et al, Astrophys. J. 494, 523 (1998). 
W. Rhode et al., Astropart. Phys. 4, 217 (1996). 

R. M. Baltrusaitis et al., Astrophys. J. 281, L9 (1984); Phys. Rev. D 31, 2192 (1985). 
K. S. Capelle, J. W. Cronin, G. Parente, and E. Zas, Astropart. Phys. 8, 321 (1998). 
Proc. of International Symposium on Extremely High Energy Cosmic Rays: Astro- 
physics and Future Observatories (Institute for Cosmic Ray Research, Tokyo, 1996). 
J. F. Ormes et al., in Proc. 25th International Cosmic Ray Conference (Durban, 1997), 
eds.: M. S. Potgieter et al., 5, 273; Y. Takahashi et al., in |4l|], p. 310. 
M. H. Salamon and F. W. Stecker, Astrophys. J. 493, 547 (1998). 
P. S. Coppi and F. A. Aharonian, Astrophys. J. 487, L9 (1997). 

R. Arnowitt and P. Nath, Phys. Rev. Lett. 69, 725 (1992); Phys. Rev. D 49, 1479 

(1994). 

D. J. H. Chung, E. W. Kolb, and A. Riotto, e-print |hep-ph/ 9805473 



V. Berezinsky, M. Kachelriess, and A. Vilenkin, Phys. Rev. Lett. 79, 4302 (1997). 

V. Berezinsky, P. Blasi, and A. Vilenkin, e-print |astro-ph/ 980327l| , to appear in 

Phys. Rev. D. 

F. Halzen, R. A. Vazquez, T. Stanev, and H. P. Vankov, Astropart. Phys. 3, 151 (1995). 
M. Boratav et al, eds., Nucl. Phys. B (Proc. Suppl.) 28B (1992). 
S. C. Corbato et al., Nucl. Phys. B (Proc. Suppl.) 28B, 36 (1992); M. Al-Seady et al., 
g, p. 191. 



m 



M. Teshima et al. Nucl. Phys. B (Proc. Suppl.) 28B, 169 (1992); M. Hayashida et al 
in 



41|, p. 205. 

R. Ghandi, C. Quigg, M. H. Reno, and I. Sarcevic, Astropart. Phys. 5, 81 (1996); e-print 
|hep-ph/9807264l . 

see, e.g., F. Halzen and E. Zas, Astrophys. J. 488, 669 (1997). 

see, e.g., E. Waxman and J. Bahcall, e-print |hep-ph/ 9807282] , submitted to Phys. Rev. 
D. 



15 



